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Abstract 

Introduction: The objective of this study was to investigate which genes are regulated by osteogenic protein-1 
(OP-1) in human articular chondrocytes using Affimetrix gene array, in order to understand the role of OP-1 in 
cartilage homeostasis. 

Methods: Chondrocytes enzymatically isolated from 12 normal ankle cartilage samples were cultured in high-density 
monolayers and either transfected with OP-1 antisense oligonucleotide in the presence of lipofectin or treated with 
recombinant OP-1 (100 ng/ml) for 48 hours followed by RNA isolation. Gene expression profiles were analyzed by 
HG-U133A gene chips from Affimetrix. A cut-off was chosen at 1.5-fold difference from controls. Selected gene array 
results were verified by real-time PCR and by in vitro measures of proteoglycan synthesis and signal transduction. 

Results: OP-1 controls cartilage homeostasis on multiple levels including regulation of genes responsible for 
chondrocyte cytoskeleton (cyclin D, Talinl, and Cyclin Ml), matrix production, and other anabolic pathways 
(transforming growth factor-beta (TGF-p)/ bone morphogenetic protein (BMP), insulin-like growth factor (IGF), 
vascular endothelial growth factor (VEGF), genes responsible for bone formation, and so on) as well as regulation 
of cytokines, neuromediators, and various catabolic pathways responsible for matrix degradation and cell death. In 
many of these cases, OP-1 modulated the expression of not only the ligands, but also their receptors, mediators of 
downstream signaling, kinases responsible for an activation of the pathways, binding proteins responsible for the 
inhibition of the pathways, and transcription factors that induce transcriptional responses. 

Conclusions: Gene array data strongly suggest a critical role of OP-1 in human cartilage homeostasis. OP-1 
regulates numerous metabolic pathways that are not only limited to its well-documented anabolic function, but 
also to its anti-catabolic activity. An understanding of OP-1 function in cartilage will provide strong justification for 
the application of OP-1 protein as a therapeutic treatment for cartilage regeneration and repair. 



Introduction 

Cartilage degeneration is one of the features of osteoar- 
thritis (OA). In order to identify cellular mechanisms 
that drive OA progression, it is necessary to understand 
the interplay between anabolic and catabolic processes 
responsible for cartilage homeostasis under physiological 
and pathophysiological states. Osteogenic protein-1 
(OP-1) or bone morphogenetic protein-7 (BMP-7) is 
one of the most potent growth factors for cartilage 
maintenance and repair identified thus far [1,2], A large 
number of in vivo and in vitro studies have shown a 
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high synthetic potency of human recombinant OP-1 
(rhOP-1; [2]). In earlier work, we found that the inhibi- 
tion of OP-1 gene expression by antisense oligonucleo- 
tides (ODNs) caused a significant decrease in aggrecan 
expression, aggrecan core protein synthesis, and proteo- 
glycan (PG) synthesis, which resulted in the depletion of 
PGs from the cartilage matrix [3]. These findings sug- 
gest that OP-1 plays a key role in maintenance of carti- 
lage integrity and homeostasis, but further work is 
needed to understand the mechanisms by which OP-1 
acts at the molecular level. 

In the current study, we used the Affymetrix Gene- 
Chip technology to monitor OP-1 regulation of 22,000 
genes from the human genome with specific emphasis 
on genes that are relevant to adult articular cartilage. 
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Those included matrix proteins, anabolic and catabolic 
gene products, as well as their intracellular regulators 
and receptors. Recently, applying the same methodology 
differential gene expression pattern in normal and 
OA cartilage tissue was identified [4]. These analyses 
revealed numerous interesting gene expression profiles, 
but per se did not allow elucidating cellular reaction pat- 
terns in response to defined extracellular stimuli. The 
goal of the current project was to evaluate the role OP- 
1 plays in regulating human articular cartilage homeos- 
tasis by using a gene array approach under conditions 
where endogenous OP-1 gene expression was inhibited 
by antisense ODNs ([3]; OP-1AS) or OP-1 signaling was 
activated and/or enhanced by rhOP-1. Key microarray 
findings were verified by real-time PCR and additional 
in vitro experiments of matrix synthesis and signal 
transduction. We found that OP-l/BMP-7 controls 
numerous metabolic pathways that are not limited to its 
direct anabolic or anti-catabolic function, but also 
related to cell growth, cell proliferation, differentiation, 
survival, apoptosis, and death. 

Materials and methods 

Materials 

Dulbecco's modified Eagle's medium (DMEM), fetal 
bovine serum (FBS), gentamicin, Ham's F-12, lipofectin, 
Opti-MEM, penicillin/streptomycin/fungizone (PSF), IX 
Platinum Quantitative PCR SuperMix-UDG and Super- 
Script III reverse transcriptase with oligo (dT) 12 _i8 were 
purchased from Invitrogen (Carlsbad, CA, USA). Phos- 
phorothioate ODN was custom synthesized by Oligos 
Etc. (Wilsonville, OR, USA). RNeasy mini kit, QIA 
shredder, RNase-free DNase kit and QuantiTect Primer 
Assay were purchased from Qiagen (Valencia, CA, 
USA). Real time polymerase chain reaction (PCR) pri- 
mers were custom synthesized by Integrated DNA 
Technologies (IDT), Coralville, IA, USA. 10,000 X SYBR 
Green 1 was purchased from Cambrex, Rockland, ME, 
USA. Recombinant human rhOP-1 was kindly provided 
by Stryker Biotech (Hopkinton, MA, USA). 

Isolation and culture of chondrocytes 

Full-thickness articular cartilage from the talus of the 
talocrural joint (ankle) from 12 human organ donors (age 
55 to 70 years old, Collins grade 0 to 1 [5]) and from the 
femur of the tibiofemoral joint (knee) from two human 
organ donors (age 67 and 73 years old, Collins grade 2) 
was obtained from the Gift of Hope Organ and Tissue 
Donor Network (Elmhurst, IL, USA) with Institutional 
Review Board approval and appropriate consent within 
24 hours of the donor's death. Knee cartilage was utilized 
for verification of the ankle cartilage results using real- 
time PCR. Chondrocytes were isolated by sequential 
digestion with pronase (2 mg/ml) for 60 minutes and 



collagenase P (0.25 mg/ml) overnight [6]. Chondrocytes 
were plated in high density monolayer culture (4 x 10 
cells/well in a six-well plate) and cultured for 24 hours in 
50% DMEM/50% Ham's F-12 supplemented with 10% 
FBS, 1% PSF, and gentamicin (50 ug/ml) for attachment 
prior to treatment with either antisense (OP-1 AS) or 
recombinant OP-1 (rhOP-1). Both treatments were 
administered for 48 hours in the absence of serum. 

Phosphorothioate ODNs 

Antisense ODNs were designed to be complementary to 
sequences in the 5'- and 3'-untranslated regions of the 
human OP-1 messenger RNA (mRNA) sequence 
(XM_030621, National Center for Biotechnology Infor- 
mation (NCBI)) as described [3]. All verification experi- 
ments with appropriate negative controls (sense and 
scrambled probes) were performed in a previous study 
[3]. For this study, the following antisense ODN was 
used: 5'-GGC-GAA-CGA-AAA-GGC-GAG-TGA-3' 
(position 237-257). 

Treatment groups 

Chondrocyte cultures were divided into three experimen- 
tal groups and treated for 48 hours as follows: 1) trans- 
fected with OP-1 AS in the presence of 10 ug/ml 
lipofectin [3]; 2) treated with 100 ng/ml of rhOP-1; and 3) 
culture control (no treatment, no serum). 

RNA Isolation 

Total cellular RNA was isolated using the RNeasy Mini 
Kit, following lysis of the cells with a Qia shredder [7] 
and included an on-column DNase digestion, according 
to the manufacturer's instructions (Qiagen). All samples 
were stored at -80°C until analyzed. 

Microarray and pathway analysis 

Gene expression profiles were analyzed by HG-U133A 
gene chips from Affimetrix (accession number: E-MTAB- 
571). At least 10 ug of RNA/per experimental group was 
required for analysis. Therefore, the RNA was pooled 
from donors in order to have sufficient RNA and to 
reduce donor-to-donor variations. Cells from all 12 
donors were treated with each experimental condition. 
The microarray data collection was in compliance with 
the Minimum Information About Microarray Experi- 
ments standard [8]. The quality of the RNA was checked 
by the Agilent Bioanalyzer (Agilent Technologies, Inc., 
Santa Clara, CA, USA), and the quality of the hybridiza- 
tion image was checked by the affyPLM model [9], To 
deal with the technical variation, each gene was measured 
by 11 different probes on the Affymetrix U133A microar- 
ray. A statistical model at the probe-level was used to 
identify the differentially expressed genes. To estimate 
the variance more efficiently with a small sample size, we 
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utilized an empirical Bayesian correction of the linear 
model [10]. Statistical significance was considered with a 
P-value of P < 0.001 and fold change larger than 1.5-fold 
between the treatment group and corresponding control. 
All the data analysis was conducted using the Bioconduc- 
tor/R package [11]. To interpret the biological signifi- 
cance of differentially expressed genes, a gene ontology 
analysis was conducted using DAVID/EASE [12]. 

Pathway analysis and classification by gene ontology 

Regulated genes (R > 1.5-fold, P < 0.001) were used as 
input for both analyses. The ingenuity pathway analysis 
system [13] was used to project genes onto known biolo- 
gical pathways (canonical pathways). The system deter- 
mines a significance value for each pathway based on an 
F-statistics that the input-genes occur randomly within 
this pathway. Grouping of genes was done by computing 
over-representation of regulated genes in gene ontology 
(GO) classes [14]. Statistical analysis consisted of 1) ana- 
lysis of differentially expressed genes under a single 
experimental condition in comparison to the correspond- 
ing control (up- or down-regulated in the presence of 
OP-1 antisense or rhOP-1); 2) analysis of differentially 
expressed genes when comparison is made between two 
treatments (OP-1 antisense and rhOP-1); and 3) gene 
ontology, when changes were analyzed within a family of 
genes according to their function (comparison was made 
between single treatment and control or between both 
treatments). Selected gene array results were verified 
experimentally in vitro or by real-time PCR. 

Validation experiments -quantitative real time PCR 

Selected gene array results were verified by real-time 
PCR. Superscript III reverse transcriptase with oligo 
(dT)i2_i8 was used to transcribe 4 ug of isolated total 
RNA into complementary DNA (cDNA) in a total 
volume of 20 ul according to the manufacturer's 

Table 1 Sequence of primers for quantitative real time PCR 

Primer Orientation Sequence 

18SrRNA 

fl-actin Forward 
Reverse 

BMP-2 

GAPDH Forward 
Reverse 

Gremlin-1 Forward 
Reverse 

IL-6 Forward 
Reverse 

IL-8 Forward 
Reverse 

LIF-1 Forward 
Reverse 



instructions (Invitrogen). Real time PCR primer sets spe- 
cific for human P-actin, GAPDH, gremlin-1, IL-6, IL-8, 
and LIF-1 (Table 1) were designed using the Primer- 
Quest program (Integrated DNA Technologies, Inc., 
Coralville, Iowa, USA). The specificity of the primers 
was verified by testing in BLAST searches [15]. Real 
time PCR primer sets specific for human 18SrRNA and 
BMP-2 were purchased from Qiagen. Real time PCR 
was performed using the Smart Cycler System (Cepheid, 
Sunnyvale, CA, USA). Each 50 ul reaction mixture con- 
tained IX Platinum Quantitative PCR SuperMix-UDG, 
0.5X Smart Cycler additive reagent (0.1 mM Tris, pH 
8.0; 0.1 mg of bovine serum albumin per ml, 75 mM 
trehalose, and 0.1% Tween 20), 0.5X SYBR Green 1 
(vendor stock 10.000X; Cambrex, Rockland, ME), 0.2 
uM each of forward and reverse primer (IDT primers) 
or IX QuantiTect primers (Qiagen primers) and 1 ul 
cDNA (WSrRNA, P-actin, BMP-2, GAPDH, gremlin-1, IL- 
6, IL-8) or 2 ul cDNA {LIF-1). Cycling parameters were: 
preheat at 60°C for 120 seconds then 95°C for 120 sec- 
onds followed by 40 three-step cycles of 95°C for 15 sec- 
onds, various annealing temperatures and times (Table 1) 
and 72°C for 30 seconds. After the last amplification 
cycle, PCR products were analyzed by melting curve ana- 
lysis in the Smart Cycler by slowly increasing the tem- 
perature to 95°C. The reactions were run in triplicate 
with appropriate controls (no cDNA template). The data 
were analyzed by using the Cepheid Smart Cycler soft- 
ware (version 2.0c) and reported as threshold cycle (C t ). 
Change in gene expression was calculated as fold change 
= 2" A(ACt) , where A(AQ) = (C t sample - C t housekeeping 
gene) - (C t control - C t housekeeping gene). 

Statistical analysis for real-time 

PCR Data are expressed as mean +/- standard deviation. 
Statistical significance was assessed by the Student i-test 
and P-values < 0.05 were considered significant. 



Annealing temp and time Accession no. 



Qiagen QuantiTect Primer Assay 
5'-TCCATCATGAAGTGTGACGTGGAC-3' 
5'-TTGATCTTCATTGTGCTGGGTGCC-3' 
Qiagen QuantiTect Primer Assay 
5'-TGGACTCCACGAGGTACTGAG-3' 
5'-CGGGAAGCTTGTGATCAATGGAA-3' 
5'-ATACCTGAAGCGAGACTGGTGCAA-3' 
5'-AACAGAAGCGGTTGATGATGGTGC-3' 
5'-GTCAATTCGTTCTGAAGAGGTGAGT-3' 
5'-CCGCAGGAGAAGATTCCAAAGATG-3' 
5'-AGACATACTCCAAACCTTTCCACCC-3' 
5'-ATTTCTGTGTTGGCGGAGTGTGGT-3' 
5'-TAAGGAGGCCTCGCAGGATGTG-3' 
5'-TAGTGGTGTACCTTGGCACCTC-3' 



62°C, 40 sec 
62°C, 40 sec 

60°C, 40 sec 
62°C, 40 sec 

64°C, 40 sec 

64°C, 40 sec 

58°C, 30 sec 

64°C, 30 sec 



[GenBank:X03205] 
[GenBank:NM_001101] 

[GenBank::NM_001200] 
[GenBank:NM_002046] 

[GenBank:NM_013372] 

[GenBank:NM_000600] 

[GenBank:NM_000584] 

[GenBank:NM_002309] 
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Results 

Microarray analysis: overview of data 

GeneChip (HG-U133A) expression data from un-stimu- 
lated, rhOP-1 and OP-1AS treated chondrocytes 
maintained in high-density monolayer culture were gen- 
erated. For the analysis of the expression data we used a 
three step analytical strategy: (I) processing of raw inten- 
sity values and normalization of profiles, (II) examina- 
tion of expression levels of gene categories that are 
relevant to articular cartilage, and (III) comparison of 
gene expression changes between the two treatments - 
OP-1AS to knockdown endogenous OP-1 expression vs. 
addition of exogenous rhOP-1. 

Analyzing the number of differentially expressed genes 
(fold changes of larger than 1.5 and corresponding 
P- values < 0.001 compared to control) after rhOP-1 or 
OP-1AS, we found that rhOP-1 modulated expression of 
4,057 genes, while OP-1AS treatment modulated expres- 
sion of only 2,618 genes respectively. More genes were 
down-regulated than up-regulated by either treatment: 
rhOP-1 down-regulated 3,365 genes vs 692 genes that 
were up-regulated; while OP- IAS down-regulated 2,364 
genes and up-regulated only 254 genes. The functional 
groups of genes modulated by lack or excess of OP-1 are 
depicted in Figure 1. RhOP-1 primarily controlled genes 
responsible for molecular function, biological processes, 
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Figure 1 Schematic representation of genes grouped accordii 

OP-1; B, genes down-regulated by treatment with recombinant OP 
regulated by OP-1 antisense treatment. 



and cellular components, while OP-1AS primarily 
affected genes controlling cellular processes and catalytic 
activity. Interestingly, either treatment up-regulated fewer 
functional groups than the number that were down-regu- 
lated (Figure 1). For example, rhOP-1 induced only five 
functional groups vs four induced by OP-1AS; while 
rhOP-1 down-regulated 19 functional groups vs 12 
down-regulated by OP-1AS. When the results were com- 
pared between the two treatments, we found that very 
few gene groups with the same function were differen- 
tially regulated by both treatments (Figure 1). Groups 
regulated by both OP-1 conditions included the genes 
responsible for cellular processes (the same number of 
genes were up-regulated by either treatment, 100 vs 101), 
development, protein binding, signal transducer activity 
and signal transduction. 

Analysis of catabolic genes: cytokines and their regulators 

Previously, we showed that OP-1 was able to counteract 
the catabolic activity of IL-ip [16,17] and other catabolic 
mediators such as fragments of cartilage matrix, fibro- 
nectin and hyaluronan [17-20]. Therefore, it was of 
interest to determine the effects of OP-1 on genes regu- 
lating pro-catabolic activity. Consistent with an anti- 
catabolic function for OP-1, a broad spectrum of genes 
with various pro-catabolic activities (cytokines and their 
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regulators, matrix degrading proteinases, apoptosis- 
related genes, neuromediators, transcription factors, and 
so on) were modulated by OP-1. Multiple cytokines and 
chemokines, in particular members of the IL-6 family, 
(Figure 2), as well as their receptors and regulators of 
their activity (Tables 2 and 3) were found to be regu- 
lated by OP-1. Interestingly, among these mediators 
only members of the IL-6 family (leukemia inhibitory 
factor (LIF), IL-11, IL-8, and IL-6) were differentially 
regulated by the two treatment conditions: rhOP-1 
down-regulated LIF expression by more than 15-fold, 
IL-11 expression by more than eight-fold, IL-8 gene by 
four-fold and IL-6 by two-fold, respectively (Figure 2A). 
Likewise, when endogenous OP-1 was inhibited by OP- 
1AS, expression of these four chemokines was elevated 
by about two-fold indicating a tight association between 
OP-1 levels and expression of members of the IL-6 
family. Verification experiments of gene array findings 



included both real-time PCR analysis and in vitro meta- 
bolic tests (Figure 2). These tests confirmed that when 
chondrocytes in high-density monolayer cultures were 
treated with rhOP-1 for 48 hours, gene expression of 
LIF, IL-6, and IL-8 was inhibited as detected by real- 
time PCR, although the magnitude of changes was dif- 
ferent from those identified by gene array (Figure 2A, 
B). In metabolic studies, we also found that OP-1 could 
overcome an inhibitory effect of IL-6 on PG synthesis in 
chondrocytes cultured in alginate beads (Figure 2C). In 
addition, our previous studies showed an ability of OP-1 
to inhibit mRNA expression of IL-1, IL-6, IL-8, and 
other cytokines in primary and immortalized chondro- 
cytes [17]. 

In analyzing the relationship between treatments to 
modulate OP-1 and the expression of genes in the IL-6 
signaling pathway, we found that OP-1 not only regu- 
lates expression of the IL-6 family of cytokines but also 
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Figure 2 Association between OP-1 and IL-6 family of chemokines. A, Effect of lack (OP-1 antisense oligo) or excess of OP-1 (treatment with 
recombinant protein, 100 ng/ml, 48 hours) on gene expression of IL-6, IL-8, IL-1 1, and LIF in chondrocytes cultured in monolayers. Graphical 
representation of gene array data. B, Real time PCR of in vitro verification experiments, where knee chondrocytes cultured in monolayers were 
treated for 48 hours with the same dose of recombinant OP-1. The graph illustrates an inhibition of LIF, IL-6, and IL-8 gene expression. C, 
verification experiments with metabolic study. Proteoglycan synthesis measured in chondrocytes cultured in alginate beads and treated for 48 
hours with 100 ng/ml IL-6 (in the presence of 150 ng/ml soluble IL-6 receptor) or the combination of IL-6 and OP-1 (100 ng/ml). Data were 
normalized to the DNA content and compared to 10% FBS control. OP-1 was able to overcome an inhibitory effect of IL-6 on PG synthesis. 
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Table 2 Changes in chemokines, cytoki 


nes, and their receptors 






Gene 


rhOP-1 vs Cntr 


OP-1 AS vs Cntr 






IUIU Lllaliyc 


IUIU Llldliytr 


MLLcbblUll I1U. 


UF 


15.86| 


2.26T 


[GenBank:NM_002309] 


IL-11 


8.694 


1.82T 


[GenBank:NM_000641] 


IL-8 


4.0 U 


1.80T 


[GenBank:NM_000584] 


IL-6 


2.09| 


1.60T 


[GenBank:NM_000600] 


IL-SRa 


2.47t 


2.40T 


[GenBank:NM_000564] 


TNF-a induced 




2.144 


[GenBank:NM_007115] 


protein 6 








TNF-a induced 


2.024 


1.60f 


LGe n Ba n k: N M_006290J 


protein-3 








TNF-RI2 


1 791 




[GenBank:NM_016639] 


TNF-R9 


1 73 


1 574, 


L^jtri 1 Da 1 1 r\.i \ 1 v 1 uu 1 ju 1 j 


TNF-R5 


1 .881 




[GenBank:NM_001 250J 


TNF-13 




1 694 


L^ci luai ii\.i\ivi uujouoj 


IL1-R1 




1.594 


[GenBank:NM_000877] 


TNF-R1 18 (osteoprotegerin) 




1.584. 


[GenBank:NM_002546] 


IL-13Rai 




1.554 


[GenBank:NM_001560] 


IL-120 


1.74| 




[GenBank:NM_002187] 


IL-1R accessory 


1.64| 




[GenBank:NM_014271] 


protein-like / 








TNF-R6 




2.084 


[GenBank:NM_000043] 



controls expression of their receptors and downstream 
intracellular mediators including signal transducers and 
activators of transcription (STATs), mitogen activated 
protein (MAP) kinases, and transcription factors. This 
suggests OP-1 inhibits IL-6 signaling at multiple levels 
(Table 3). Among other genes that either regulate cyto- 
kine activity or mediate their signaling, the most affected 
by OP-1 were the receptors for IL-ip and tumor necro- 
sis factor alpha (TNF-a) (see Table 2) as well as TNF-a 
inducible protein. Although under the experimental 
conditions expression of TNF-a and IL-ip genes was 
not influenced by OP-1, previous studies showed that 
injection of OP-1 into nucleus pulposus inhibited 



production of autocrine TNF-a and IL-ip elevated in 
response to injurious compression of the intervertebral 
discs [21] proving an association between OP-1 and sig- 
naling pathways of the above mentioned cytokines. In 
addition, several other studies have provided evidence of 
an ability of OP-1 to regulate either IL-ip induced 
responses or IL-ip downstream signaling [16-18,22,23]. 

Analysis of catabolic genes. Neuromediators 

Previous studies have provided evidence that OP-1 may 
regulate mediators of pain-related behavior and their 
activation in response to injurious compression of the 
intervertebral disc and acute cartilage trauma [24-26]. 



Table 3 Changes in the mediators of IL-6 signaling pathway 



rhOP-1 vs Cntr 
fold change 



OP-1 AS vs Cntr 
fold change 



Accession no. 



Genes from IL-6 signaling pathway 

ELK-1 

IL-6 

IL-6R 

IL-6 signal transducer (oncostatin M receptor) 

SX477 

NFkBIoc 

Protein inhibitor of activated STAT3 
STAT6 

MAP 3 kinase 7 
MARK 14 
MAPK1 



2.094 
1.81 4 

2.424 



1.60| 

1.634 

58t 
844 
534 
674 
524 



[GenBank:NM_005229] 
[GenBank:NM_000600] 
[GenBank:NM_000565] 
[GenBank:NM_002184] 
[GenBank:NM_007315] 
[GenBank:NM_020529] 
[GenBank:NM_006099] 
[GenBank:NM_003153] 
[GenBank:NM_003188] 
[GenBank:L35253] 
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We also reported that injection of OP-1 into nucleus 
pulposus down-regulated substance P expression [26], 
bradykinin and bradykinin inducible receptor pi [26]. 
Therefore, it was of interest to examine expression of 
neuromediators and their receptors in the present array 
study. After stimulation for 48 hours with rhOP-1, 
expression of the receptors of bradykinin and substance 
P was down-regulated (Table 4). Both receptors of bra- 
dykinin, constitutively expressed P2 and inducible pi, 
were down-regulated by the treatment with OP-1. 
Expression of the pi receptor was differentially regu- 
lated under conditions of excess and lack of OP-1, that 
is, treatment with rhOP-1 inhibited gene expression of 
this receptor by 1.85-fold, while its expression was up- 
regulated by 1.59-fold when endogenous OP-1 expres- 
sion was inhibited by antisense oligonucleotides. These 
results are consistent with previous data on the protein 
level in an in vivo model of disc herniation, where injec- 
tion of OP-1 into the nucleus pulposus completely abol- 
ished bradykinin receptor P 1 [26] . Although by gene array 
we did not identify significant changes in the expression of 
bradykinin and substance P at the time point tested here, 
we found changes in substance P receptor and its precur- 
sor. We also found that OP-1 inhibited expression of 
nerve growth factor-P by almost two-fold. 

Analysis of catabolic genes: Transcription factors 

Besides cytokines and their receptors, OP-1 also affected 
gene expression of transcription factors that regulate 
cytokine signaling. Previously, in normal primary and 
immortalized chondrocytes, we found that OP-1 inhibits 
activation of the nuclear factor kappa-light-chain-enhan- 
cer of activated B cells (NF-kB) and activator protein- 1 
(AP-1) transcription factors [17]. Here, expression of a 
large set of transcription factors was found to be modu- 
lated by OP-1 (Table 5). In addition to common factors 
such as NF-kB, STAT1 and STAT6, gene array also dis- 
covered factors that repress IL-2 expression, p38 inter- 
acting protein, Runxl, and others. The majority of these 
transcription factors regulate directly or indirectly (as 
p38 interacting protein) transcriptional responses 
induced by various pro-inflammatory mediators (IL-ip, 
IL-6, matrix fragments). Others, like Runxl, are involved 

Table 4 Changes in neuromediators and their receptors 



in the process of chondrogenesis. To further demon- 
strate the effect of OP-1 on activation of transcription 
factors, we treated cultured cells and found that OP-1 
was able to at least partially inhibit activation of NF-kB 
in primary chondrocytes pre-treated with IL-ip or acti- 
vation of Stat-1 in chondrocytes treated with IL-6 and 
IL-6 soluble receptor (data not shown). 

Analysis of catabolic genes: Matrix degrading proteases, 
cathepsins, and apoptosis-related genes 

Among other catabolic genes influenced by OP-1 were 
the matrix metalloproteinases (MMPs), cathepsins, and 
a number of proteases with various modes of action 
(Table 6). Thus, expression of membrane type-1 MMP 
{MMP-14) was inhibited by rhOP-1 by 1.6-fold (P < 
0.001) along with tissue inhibitor of metalloproteinases 
(TIMP)-3 (2.06-fold, P < 0.001). At the same time, 
expression of MMP-2 (gelatinase A), which is activated 
by MMP-14 [24], was not affected by rhOP-1, but was 
down-regulated by OP-1AS (2.31-fold, (P < 0.001) as 
well as was MMP-9 (gelatinase B) (1.5-fold). Interest- 
ingly, the same positive correlation was found between 
the levels of OP-1 and expression of another TIMP, 
TIMP-4, which was decreased by 1.7-fold in the OP- 
1AS group confirming its association with MMP-2 [25]. 
Parallel changes were observed in other types of pro- 
teases, such as a disintegrin and metalloproteinases 
(ADAM)-9, 10, and 28. Their gene expression was down- 
regulated under OP-1AS from 2.34 to 1.75-fold. Treat- 
ment of chondrocytes with rhOP-1 inhibited expression 
of ADAM-15,-19, as well as urokinase type plasminogen 
activator, its receptor, and transglutamianse-2. There 
were also some proteinases that were up-regulated by 
rhOP-1: ADAM-TS7, ADAM-TS12, and tissue specific 
plasminogen activator suggesting that perhaps these pro- 
teins are involved in anabolic/remodeling processes. 

Among the proteases that were also regulated by OP-1 
were cathepsins B, C, and S. So far, these lysosomal 
cysteine proteases have been less studied in cartilage, 
though cathepsin C appears to be a central coordinator 
for activation of many serine proteases in immune/ 
inflammatory cells [29], while cathepsin B was thought 
to play an important role in the development of 



rhOP-1 vs Cntr OP-1 AS vs Cntr 

fold change fold change Accession no. 

Bradykinin Rpi 1 .851 1 .59t [GenBank:NM_000710] 

Bradykinin RP2 1.68 J, [GenBank:NM_000623] 

Tachykinin R1 1.641 [GenBank:NM_001058] 

Nerve growth factor-p 1.931 [GenBank:NM_002506] 

Tachykinin 1 precursor 2.261 
(Substance K, Substance P) 
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Table 5 Changes in transcription factors 





rhOP-1 vs Cntr 


rhOP-1AS vs Cntr 






fold change 


fold change 


Accession no. 


Transcription factor 8 (represses IL-2 expression) 


3.284 


2 974 


[GenBank:NM_030751] 


NF-K.B2 


2 77 1 




[GenBank:NM_002502] 


5TAT1 


2421 




[GenBank:NM_00731 5] 


Transcription factor AP-2a 


2.074 


1.524 


[GenBank:NM_003220] 


Suppression of tumorigenicity 


2.044 


2.204 


[GenBank:NM_01 3437] 


Runx 1 


1.89 J. 


1.644 


[GenBankNM 001754] 


NFnBia 


1.86 J, 


1.58T 


[GenBank:NM_020529] 


NFYfi 


1.681 


1.754 


[GenBank:NM_006166] 


Activating transcription factor 7 


1.66| 




[GenBank:NM_006856] 


MADS box transcription enhancer factor 2-d 


1.654 




[GenBank:NM_005920] 


Upstream transcription factor 2, c-fos interacting 


1.584 




[GenBank:NM_003367] 


Transcription factor (p38 interacting protein) 


1.574 




[GenBank:NM_017569] 


MADS box transcription enhancer factor 2-C 


1.60T 


1.944 


[GenBank:NM_002397] 


Protein inhibitor of activated STAT3 




1.844 


[GenBank:NM_006099] 


Ubiquitin-like 1 (sentrin) 




1.604 


[GenBank:NM_003352] 


STAT6 




1.534 


[GenBank:NM_003153] 


Table 6 Changes in proteases and their inhibitors 




rhOP-1 vs Cntr 


rhOP-1AS vs Cntr 






fold change 


fold change 


Accession no. 


Bcl-2 


2.454 




[GenBank:NM_001191] 


Caspase 4, apoptosis-related cysteine protease 


2.114 


1.594 


[GenBank:NM_001225] 


Programmed cell death 8 (apoptosis-inducing factor) 


1.704 




[GenBank:NM_004208] 


Calpain 9 


1.554 




[GenBank:NM_006615] 


Caspase 6 




2.184 


[GenBank:NM_001226] 


Caspase 8 




1.824 


[GenBank:NM_001228] 


Caspase 2 




1.50| 


[GenBank:NM_001224] 


MMPs and inhibitors 








TIMP-3 


2.064 




[GenBank:NM_000362] 


MMP-14 


1.554 




[GenBank:NM_004995] 


MMP-2 




2.314 


[GenBank:NM_004530] 


TIMP-4 




1.694 


[GenBank:NM_003256] 


MMP-9 




1.504 


[GenBank:NM_004994] 


ADAM and ADAMTS 








ADAM-19 


1.834 




[GenBank:NM_023038] 


ADAM-15 


1.514 




[GenBank:NM_003815] 


ADAMTS- 12 


1.88T 


2.034 


[GenBank:NM_030955] 


ADAMTS-7 


1.58| 




[GenBank:NM_014272] 


ADAM-10 




2.344 


[GenBank:NM_001110] 


ADAM-28 




1.634 


[GenBank:NM_014265] 


ADAM-9 




1.614 


[GenBank:NM_003816] 


ADAM-7 




1.564 


[GenBank:NM_003817] 


Cathepsins 








Cathepsin B 




2.144 


[GenBank:NM_001908] 


Cathepsin C 




1.754 


[GenBank:NM_001814] 


Cathepsin S 




1.754 


[GenBank:NM_004079] 


Other proteases 








Transglutaminase 2 


2.104 




[GenBank:NM_004613] 


Plasminogen activator-urokinase 


1.574 




[GenBank:NM_002658] 


Tissue Plasminogen Activator 


1.564 




[GenBank:NM_000930] 



Chubinskaya ef al. Arthritis Research & Therapy 201 1, 13:R55 
http://arthritis-research.eom/content/13/2/R55 



Page 9 of 14 



osteoarthritis [30]. Expression of all three cathepsin 
genes was down-regulated under OP- IAS. 

A previous study on acute impact injury in vivo [31] 
strongly suggested an anti-apoptotic effect of OP-1 in 
post-traumatic OA. Therefore, we expected that OP-1 
may control genes involved in apoptosis-related pro- 
cesses. We found that rhOP-1 inhibited program cell 
death 8 gene (apoptosis-induced factor), Bcl-2 gene and 
the calpain-9 gene (Table 6). However, the key caspases 
that trigger and promote cell death by apoptosis were 
not affected. During the absence of OP-1 (antisense 
treatment), expression of caspases 8, 9, and 6 were 
inhibited and only caspase 2 was elevated (Table 6). The 
reason for a down-regulation of the apoptosis-related 
genes under conditions where OP-1 is lacking is not 
clear, but may be a response to help avoid cell death. 

Analysis of anabolic genes: transforming growth factor- 
beta (TGF-p)/BMP family, their receptors and regulators of 
signaling 

Affimetrix analysis identified a very interesting effect of 
OP-1 on members of the BMP/TGF-|3 family (Table 7). 



Treatment with rhOP-1 down-regulated expression of 
growth differentiation factor (GDF)-15, BMP-2, and Acti- 
vin A, and BMP-2 inducible kinase, while inhibition of OP- 
1 expression up-regulated GDF-15 and Activin A. Down- 
regulation of BMP-2 expression in chondrocytes treated 
with rhOP-1 was confirmed by real-time PCR (Figure 3). 
Antisense reduction of OP-1 levels resulted in down-regu- 
lation of GDF-10 and TGF-a expression (Table 7). Further- 
more, a correlation was also found between OP-1 and the 
mediators of its downstream signaling, where OP-1AS 
treatment inhibited expression of transcription factors, Id 
proteins 2 to 4 (Table 7), binding protein Gremlin (Figure 
2), and MAD genes. Changes in Id genes correlated with 
the earlier findings from our laboratory, which demon- 
strated that the treatment of chondrocytes with rhOP-1 led 
to the elevation of Idl, Id2, and Id3 genes and proteins 
[32]. Contrary to changes in the Gremlin gene, which 
showed a positive correlation with OP-1 levels, expression 
of Follistatin binding protein was inhibited by more than 
two-fold in chondrocytes treated with rhOP-1. 

In addition, OP-1 modulated expression of the TGF- 
P/BMP receptors. With the exception of Activin-a RIB, 



Table 7 Changes in the expression of TGF-p7BMP family related genes, their receptors, and signaling regulators 



rhOP-1 vs Cntr rhOP-IAS vs Cntr 

fold change fold change Accession no. 



GDF-15 


3.044 


2.03T 


[GenBank:NM_004864] 


BMP-2 


2.67| 




[GenBank:NM_001200] 


lnhibin-pa (activin A) 


2.324 


2.15| 


[GenBank:NM_002192] 


BMP-2 inducible kinase 


1.614 




[GenBank:NM_017593] 


Parathyroid hormone-like hormone 


1.604 


2.1 7f 


[GenBank:NM_002820] 


ID2 




2.324 


[GenBank:NM_002166] 


Notch 4 




2.324 


[GenBank:NM_004557] 


MAD-6 




2.054 


[GenBank:NM_005585] 


Gremlin 


1.88T 


1.944 


[GenBank:NM_013372] 


GDF-10 




1.864 


[GenBank:NM_004962] 


ID4 


1.82T 


1.844 


[GenBank:NM_001546] 


ID3 




1.734 


[GenBank:NM_002167] 


MAD interacting protein 




1.694 


[GenBank:NM_004799] 


MAD-4 




1.674 


[GenBank:NM_005359] 


Notch 1 




1.654 


[GenBank:NM_017617] 


MAD-7 




1.624 


[GenBank:NM_005904] 


TGF-a 




1.544 




Receptors 








Frizzled homolog 10 (Drosophila) 


1.574 




[GenBank:NM_007197] 


Activin-a Rl 


1.534 




[GenBank:NM_001105] 


Activin A-R1IB 




2.424 


[GenBank:NM_001106] 


BMPR1A 




1.834 


[GenBank:NM_004329] 


TGF-f)RI 




1.514 


[GenBank:NM_004612] 


TGF-pRIII 




1.504 


[GenBank:NM_003243] 


TGF-P R2 


1.58T 




[GenBank:NM_003242] 


Activin A-RIB 




1.53f 


[GenBank:NM_004302] 


Bone formation 








Osteomodulin 


1.78t 


2.564 


[GenBank:NM_005014] 
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OP-1 treated 





P<0.001 






















3.98-fold 















GAPDH 18SrRNA BMP-2 

Gene names 

Figure 3 Effect of OP-1 on BMP-2 gene expression. Real time 
PCR of in vitro verification experiments, where knee chondrocytes 
cultured in monolayers were treated for 48 hours with 100 ng/ml 
recombinant OP-1. The graph illustrates an inhibition of BMP-2 
mRNA expression. 



which was inhibited by rhOP-1 and elevated under the 
lack of OP-1, expression of other receptors, Activin-a 
RUB, BMPRIA, TGF-B RI, II, and III correlated posi- 
tively with OP-1 expression (Table 7). 

Analysis of anabolic genes: other growth factors 

Previously we showed that rhOP-1 stimulated expres- 
sion of insulin-like growth factor (IGF)-l and IGF-1 
receptor genes [17], while inhibition of OP-1 gene 
expression by OP- IAS down-regulated mRNA expres- 
sion of these genes. We have also documented a syner- 
gistic effect of OP-1 on IGF-1 induced responses in 
normal and OA chondrocytes [33,34]. Here, we con- 
firmed an association between OP-1 and IGF-1 path- 
ways by documenting a 1.73-fold decrease in IGF-1 
receptor expression and a decrease in two IGF-1 binding 
proteins-5 and 7 (1.9- and 1.5-fold respectively) under 
OP-1AS. Furthermore, other genes within the IGF-1 sig- 
naling pathway were regulated by OP-1. Among them 
were PIK3R1, PRKAR2B, MAP2K2, PDE3B, and SOCS3 
(Table 8). 

Modulation of OP-1 levels affected mRNA expression 
of growth factors and some of their receptors that belong 
to various families, such as Nerve Growth Factor-p, Vas- 
cular Endothelial Growth Factor, Endothelial Cell Growth 
Factor 1 (platelet-derived), Capillary Morphogenesis Pro- 
tein-1, and Fibroblast Growth Factor (FGF)-7. Their 
expression was inhibited by rhOP-1 from 1.93- to 1.5- 
fold. Contrary, the expression of the FGF-R2 and 3 recep- 
tors, and a and P receptors of Platelet-Derived Growth 
Factor was stimulated by rhOP-1 Table 8). 

Matrix proteins and their receptors 

Cartilage-specific matrix genes were found to be 
modulated by rhOP-1 treatment. Expression of the 



Table 8 Association between OP-1 and other growth 
factors including igf-1, insulin, and tyrosine-kinase 
signaling 





rhOP-1 vs 


OP-1 AS vs 




Cntr 


Cntr 




told 


fold change 




change 




IGF-BPi 


2 171 




Mprvp nrnwth fsrfnr-Pt 


1 .93J, 




VEGF-B 


1 .62 J, 


1 .504 


Endothelial cell growth factor 1 


1.561 




(platelet-derived) 






VEGF 


1.524 




Capillary morphogenesis protein 1 


1.52| 




FGF-7 




2.874 


FGF-R2 


1.69t 


2.834 


IGF-BP5 




1.904 


FGF-H3 


1.87t 


1.804 


iCC 7 D 

ihr-ln 




1 .734 


PDGF-Ra 


1.62| 


1.704 


PDGF-Rp 




1.684 


IGF-BP7 




1.584 


IRS2 (insulin receptor substrate 2) 


2.104 


1.704 


DPYSL2 (dihydropyrimidinase-like 2) 


1.604 


1.604 


MET (hepatocyte growth factor receptor) 


1.704 


1.604 


SPRY2: sprouty homolog 2 (Drosophila) 


1.604 


1.604 


S0RBS1: sorbin and SH3 domain containing 
1 


1.704 


1.504 


PIK3R1 (Phosphoinositide-3-kinase, 


1.724 




regulatory subunit 1) 






MAP2K2 (mitogen-activated protein kinase 


1.614 




kinase 2) 






PDE3B (phosphodiesterase 3B, cGMP- 


2.004 




inhibited) 






SOG53 (suppressor of cytokine signaling 3) 


1.794 





collagen type IX-a3 chain and cartilage oligomeric 
protein (COMP) was up-regulated by about 1.5-fold in 
chondrocytes treated with rhOP-1 (Table 9). Among 
proteoglycans, versican was affected the most (by 
about three-fold down-regulation by OP- IAS) and syn- 
decan was differentially regulated under both rhOP-1 
and OP- IAS treatments. There were a number of 
other matrix genes regulated by OP-1: bone sialopro- 
tein, osteonectin, cadherins, chondroitin sulfate PG4 
and dermatan sulfate PG3 (Table 9). As expected, 
there was a positive correlation between OP-1 and 
CD44 expression. Inhibition of OP-1 expression 
resulted in 2.34-fold reduction in CD44 expression. 
However, contrary to previously published data [35], 
rhOP-1 inhibited hyaluronan synthase 2 expression. 
A number of basement membrane proteins were 
modulated by OP-1: al,2,3, and five chains of collagen 
type IV, laminin, versican among others. Gene expres- 
sion of bamacan and laminin was inhibited by rhOP-1 
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Table 9 Changes in the expression of matrix proteins, their receptors, and integrins 

rhOP-1 vs Cntr OP-1AS vs Cntr 

fold change fold change Accession no. 



Matrix proteins 




Collagen IV-ct3 


1.82| 


Laminin-pi 


1 .654. 


Chondroitin sulfate PG6 (bamacan) 


1.60| 


Versican (chondroitin sulfate PG2) 




Collagen l-ctl 




Collagen XIV-ct1 




Collagen l-ct2 




Cartilage associated protein 




Cadherin 1 1 (OB-cadherin (osteoblast)) 




Collagen XVI-or.1 




Dermatan Sulfate PG3 


1.57t 


Collagen V-al 




Bone sialoprotein 




Collagen Vlll-a2 




Collagen Vl-a1 




Collagen IV-a1 




Collagen IV-a2 




Syndecan 


1.58| 


Collagen V-a2 




Osteonectin 




Cadherin 19 




/~^ll— I \/ 

Collagen IX-a3 


1 .59^ 


Cadherin 


1.54| 


COMP 


1 52\ 


Collagen IV-oc5 




Chondroitin sulfate PG4 




Matrix protein receptors 




HAS2 


1.78| 


CD44 




Integrins 




htegrin-a5 


1.77| 


Integrin-p4 


1.644 


htegrin-a6 




Integrin-p-llke 1 




Integrin-p3 




htegrin-aE 





and stimulated under OP-1AS. OP-1 also modulated 
expression of collagens that are not cartilage-specific, 
such as collagen type I, IV, V, VI, VIII, XIV, and XVI. 
Their expression was inhibited under the OP-1AS 
treatment (Table 9). The greatest decrease in mRNA 
expression was found for al and cx.2 chains of type I 
collagen (more than 2.6-fold). 

Discussion 

To the best of our knowledge, this is the first report that 
uses microarray analysis to provide a comprehensive 
understanding of the role OP-1 plays in overall cartilage 



1.544 [GenBank:NM_000091] 

3.384 [GenBank:NM_002291] 

1.52| [GenBank:NM_005445] 

2.984 [GenBank:NM_004385] 

2.634 [GenBank:NM_000088] 

2.594 [GenBank:NM_021110] 

2.574 [GenBank:NM_000089] 

2.13| [GenBank:NM_006371] 

2.03| [GenBank:NM_001797] 

2.00| [GenBank:NM_001856] 

1.954 [GenBank:NM_004950] 

1.894 [GenBank:NM_000093] 

1.834 [GenBank:NM_004967] 

1.734 [GenBank:NM_005202] 

1.704 [GenBank:NM_001848] 

1.704 [GenBank:NM_001845] 

1.694 [GenBank:NM_001846] 

1.674 [GenBank:NM_002997] 

1.654 [GenBank:NM_000393] 

1.514 [GenBank:NM_003118] 

1.504 [GenBank:NM_021153] 
[GenBank:NM_001853] 
[GenBank:NM_001408] 
[GenBank:NM_000095] 

1.89t [GenBank:NM_000495] 

1.64T [GenBank:NM_001897] 

[GenBank:NM_005328] 
2.344 [GenBank:NM_000610] 

[GenBank:NM_002205] 
[GenBank:AF011375] 
2.344 [GenBank:NM_000210] 
2.074 [GenBank:NM_004791] 
1.724 [GenBank:NM_000212] 
1.704 [GenBank:NM_002208] 



homeostasis. We found that OP-1 controls cartilage 
homeostasis on multiple levels including regulation of 
genes responsible for chondrocyte cytoskeleton {cyclin D, 
Talinl, and Cyclin Ml, for example, and confirmed in 
[36]), matrix production and other anabolic pathways, as 
well as regulation of cytokines and various catabolic path- 
ways responsible for matrix degradation and cell death. 
Importantly, in many of these cases, OP-1 modulated the 
expression of not only the ligands, but also their recep- 
tors, mediators of downstream signaling, kinases respon- 
sible for an activation of the pathways and transcription 
factors that induce transcriptional responses. 
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Due to high variability among human samples, only a 
few studies have utilized microarray analysis to test the 
entire human genome in primary adult articular chon- 
drocytes [4,37-39], and only one of Saas et al. [4] 
addressed in part the effect of BMP-7/OP-1. These ana- 
lyses used the tissue from one or a maximum of two 
donor cartilage samples. In the present study, normal 
(grade 0) articular cartilage was collected from 12 
donors within a similar age range. One of the limitations 
of the study is that we examined gene expression pro- 
files only at one time point, after 48-hours of culture. 
Therefore, changes in early-response genes and late- 
response genes might have been missed. This could 
explain some results, as for example, the lack of changes 
in the expression of major cartilage matrix proteins. 
However, such an approach gave us a breath of the 
overall effects of OP-1 on cartilage homeostasis. 

Due to the abundance of the results, we will discuss 
only the most relevant and those that could be 
explained by the current knowledge of the field. Perhaps 
most important was the finding that OP-1 is a unique 
growth factor in its capacity to display simultaneously 
pro-anabolic and anti-catabolic activities. It was pre- 
viously shown that OP-1 stimulated expression and 
synthesis of collagen type II, aggrecan, hyaluronan, and 
CD44 [1,2,20,40] as well as IGF-1, IGF-1 receptor, and 
responses to IGF-1 [17]. In the current studies, we used 
high-density monolayers while in previous work explants 
or alginate beads were used with different media condi- 
tions (no serum vs serum or ITS-media). The finding 
that the microarray results shown here were comparable 
to the previous results suggest that the pro-anabolic 
effects of OP-1 in human articular chondrocytes are 
persistent. With regard to the anti-catabolic activity, the 
ability of OP-1 to counteract various pro-inflammatory/ 
catabolic responses or directly inhibit expression of cata- 
bolic mediators was previously shown in primary chon- 
drocyte cultures or in animal models of post-traumatic 
osteoarthritis or disc degeneration [17-19,24,31]. In this 
study, we found that OP-1 inhibits expression of IL-6 
and members of the IL-6 family of chemokines as well 
as their receptors and signaling mediators. Further- 
more, the tight association between these two classes 
of mediators (OP-1 and IL-6) was documented under 
both experimental conditions (plus or minus OP-1). 
Based on our new data on the role of IL-6 in acute 
post-traumatic responses [41], it is possible that OP-1 
was able to protect cartilage from degenerative 
changes caused by acute trauma [31] not only due to 
its direct effect on matrix synthesis, but also because 
of its ability to inhibit IL-6, TNF-ot, and the catabolic 
pathways induced by the fragments of the extracellular 
matrix: fibronectin [19], hyaluronan [20], and collagen 
telopeptides [42]. 



Another important effect of OP-1 may be an ability to 
inhibit expression of neuromediators and their receptors. 
Previously, an anti-pain effect of OP-1 was documented 
in the rat models of herniated disc or disc degeneration 
induced by injurious compression. In these studies, OP-1 
injections reduced hyperalgesia and inhibited elevation of 
IL-1, TNF-a, substance P, bradykinin and their receptors 
in various disc tissues including spinal cord and dorsal 
ganglion [21,24,26]. In chondrocytes, it is the first report 
that indicates a connection between OP-1 and various 
neuromediators, though substance P and its NK-1 recep- 
tor were already identified in cartilage in the model of 
mechanical stress [43]. Very recently, our findings 
received another proof in phase I placebo-controlled clin- 
ical studies on OP-1 treatment for osteoarthritis patients 
[44], in which a single injection of OP-1 reduced pain 
even after six months of treatment. 

Interesting results were found with regard to the abil- 
ity of OP-1 to regulate the TGF-p7BMP signaling path- 
way. The levels of OP-1 expression positively correlated 
with the expression of activin-like kinase (ALK)-3 or 
BMP-RIA, transcription factors Id proteins 2 to 4, and a 
binding protein Gremlin, indicating that this could be a 
primary route for OP-1 signaling. We also found that 
another binding protein, Follistatin, exhibited a negative 
correlation with the levels of OP-1. Thus, our results 
suggest a differential regulation of these two binding 
proteins by OP-1, which could mean that Gremlin and 
Follistatin perform distinct functions in the mediating 
BMP responses or they are involved in different stages 
of signaling. This point of view is supported by studies 
of Tardif et al. [45] who reported their differential 
expression and spatial distribution in the dog OA 
model. One of the most surprising results was the find- 
ing that OP-1 inhibits expression of another member of 
the BMP family, BMP-2, which shares the same signal- 
ing machinery and in many cases exhibits similar ana- 
bolic activities [23,46]. This result was confirmed by 
real-time PCR and definitely warrants further studies in 
understanding the responses induced by homologous, 
yet very different members of the same family [16]. 

Finally, another unexpected result was the inhibition of 
TIMP-3 expression by rhOP-1; while previously, TGF-|3 
has been shown to induce this inhibitor [47]. The differ- 
ences in the results could be attributed to distinct culture 
conditions (primary chondrocytes in our case and pas- 
saged chondrocytes as in Qureshi et al. [47], assessment 
at various time points, or distinct signaling mechanisms 
between TGF-p* and OP-1 that are responsible for induc- 
tion of TIMP-3 expression. However, on the protein level 
it has already been reported that OP-1 inhibits TIMP-3 
protein [48]. Furthermore, changes in TIMP-3 were par- 
allel to changes in certain genes responsible for apoptosis, 
which supports the notion that in cancer cells TIMP-3 
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may promote cell death by apoptosis [28]. On the other 
hand, TIMP-3 has been shown to inhibit aggrecanase- 
mediated release of glycosaminoglycans in bovine nasal 
cartilage [49]. At this point, the role of TIMP-3 in 
human chondrocytes and its regulation by various media- 
tors remains to be investigated. 

Conclusions 

This analysis of gene array data strongly suggests a criti- 
cal role of OP-1 in human cartilage homeostasis. OP-1 
regulates numerous metabolic pathways that are not 
only limited to its anabolic function, but also to its anti- 
catabolic activity. Understanding of OP-1 function in 
cartilage will provide strong justification for the applica- 
tion of OP-1 protein as therapeutic treatment for carti- 
lage regeneration and repair. 
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